Possible Suppression of Neutron EDM 
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ABSTRACT 

Employing generalized Schiff's transformation on electric dipole moments 
(EDM) in quantum field theory, we show that the chromoelectric EDM lagrangian 
density is transformed into the electric EDM term with a new coefficient. Under 
the new constraint on the EDM operators, the neutron EDM can be described by 
a unique combination of electric EDM df and chromoelectric EDM df of quarks. 
If the special relation of df = 7^~df holds, then the neutron EDM is suppressed 
significantly. 
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1. Introduction 

Recent measurements of electric dipole moments (EDM) in neutrons and 
atoms have presented tight upper bounds [1-5] which are rather close to 
theoretical predictions of the EDM by supersymmetric theory calculations 
[6-10]. Even though there are some free parameters in the calculations, 
it is expected that the finite EDM values may well be observed soon. In 
particular, the recent measurement of the muon g — 2 [23] should give severe 
constraints on the theoretical EDM values. 

Recently, there have been many theoretical investigations of the neutron 
EDM [11-18]. In particular, Pospelov and Ritz [19] suggest that the contri- 
butions to the neutron EDM which come from the chromoelectric interac- 
tions may well be comparable to the electric EDM estimation. Therefore, 
careful studies of the EDM operator arising from the chromoelectric fields 
must be quite important for the reliable estimation of the neutron EDM. 

In this Letter, we study a constraint on the chromoelectric EDM operators 
due to the unitary transformation which is analogous to Schiff's theorem 

in nonrelativistic quantum mechanics. Here, we show that the chromoelec- 

df - 

trie EDM lagrangian density —i—if)a flu 'y 5 t a ij)G tJ '' / ' a is transformed into a 

6 f ~ — 

new term i- — dfipa^j^F^ which is just the same as the electric EDM 



lagrangian density. Thus, the coefficient of the EDM lagrangian density is 



lations should be described with the combinations of dt —dt . This 



should present some constraints which should be satisfied by any calcula- 
tions of the neutron EDM. 

2. Schiff's theorem in QED 

In 1963, Schiff proved [20] that the effect of the EDM interaction cannot be 
measured in the nonrelativistic quantum mechanics with electromagnetic 
interactions. This can be easily seen since the hamiltonian with the EDM 
interaction 
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Therefore, any neutron EDM calcu- 




H=£- + eA (r) - d E 



(1) 
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can be rewritten with the replacement of r — > r as 

e 

H = ^+eA (r)+d.VA (r) = J^-+ e A (\r+*\)+O(d 2 ) = ^- +e A (r)+O(d 2 ) 

where E(r) = — VA (r), and we have made use of the fact that the momen- 

d 

turn p does not change under the replacement of r — > r . This means 

e 

that the spectrum of the original hamiltonian with the EDM term becomes 
just the same as that of the hamiltonian without the EDM term since the 
0(d 2 ) term is negligibly small. This is Schiff's theorem, and the EDM 
effect cannot be observed in the nonrelativistic quantum mechanics with 
electromagnetic interactions. 

Next, we treat the field theory version of Schiff's theorem and consider the 
unitary transformation in QED [21]. The lagrangian density of QED with 
the electromagnetic EDM terms can be written as 

C = i;i(d, + ieA^ - m ijij - \f^ - i^a^F"* (3) 
where F^ u denotes the field tensor given as 

F, u = d,A v - d v Ar (4) 
Now, we consider the following unitary transformation with = id^, 

V>' = exp ^^75Pm7^ ^- (5) 

Under this transformation, the total lagrangian density becomes, up to the 
order of df 

C! = ^(^ + ze^) 7 ^-mo^-^F^F^-^^75(p^-e^)^- (6) 

This is the new lagrangian density in QED and states that the new EDM 
term does not couple to the external electric field while the original EDM 
term has the coupling with the electric field E in the first order as seen in 
eq.(3). This is just the same as Schiff's statement that the EDM interaction 
does not have the first order perturbation energy with the external electric 
field any more. Here, it should be noted that, even though eq.(6) contains 
the A^ term, the hamiltonian constructed from eq.(6) contains only the 
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three dimensional vector potential A ■ p term which does not contribute 
to the EDM. This is wellknown in the EDM evaluation in atomic physics 
[1,21,25]. 

Therefore, the EDM of a composite system can be generated in the second 
order perturbation as 

, _ < N\E i e t z t \N* >< N*\H edm \N > 

«n — * {() 

where the sum should be taken over excited states in the composite system 
and H edm corresponds to the EDM lagrangian term in eq.(6). N and N* 
denote the ground and excited states of the composite system. 

At this point, we make a comment on the gauge sysmmetry violation of the 
unitary transformation of eq.(5). This transformation leads to the violation 
of the gauge symmetry, but this is justified as long as one is interested in 
evaluating the low energy state property of the system up to the order of g. 
This is the same as Weinberg-Salam treatment of spontaneous symmetry 
breaking where the gauge symmetry is broken when evaluating the ground 
state property of the system. Even though Weinberg-Salam's model is not 
exact, their treatment is economical. Clearly, one can obtain the same 
result with the spontaneous symmetry breaking procedure as the result 
that is solved exactly for the low energy property. But obviously the exact 
calculation of the Weinberg-Salam's model is far from economical. 



3. Unitary transformation in QCD plus QED 

Recently, there have been several papers which treat the chromoelectric 
EDM terms [15-19]. The neutron EDM which comes from the chromo- 
electric EDM of quarks seems to depend on the model calculations. The 
calculations with QCD sum rules in ref. [16] suggest a suppression of the 
neutron EDM from the chromoelectric EDM of quarks while the chiral loop 
estimate gives a sizable contribution to the neutron EDM [18]. Further, 
Pespelov and Ritz [19] show that the chromoelectric EDM term gives a 
significant contribution to the neutron EDM when they employ QCD sum 
rules. 
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In this section, we show that the dominant contributions from the chromo- 
electric EDM can be transformed into the identical shape to the electric 
EDM term with some modified coefficients. 

We start from QCD plus QED lagrangian density with electric and chro- 
moelectric EDM terms where the charge of the / flavor quark is denoted 
by e f , 

C = Md, + tg s A;t a + iefAjW - m ^ - \f, v F^ - \G^G^ a 

-^CV75*W^ - i^o^F^ (8) 

where 

G% = d,A a v - d v A% - g s C ahc A\A c v . (9) 
Now, we consider the following unitary transformation 

' = 6XP (^ 75 °° M " 9sA V a)r ) (10) 
Under this transformation, we obtain the following new lagrangian density, 

£ = Md^ + ig s Ap a + ie f A^ - m ^ - \f^ - \c%G^ a 

-\ (d f - ^-df^J ^Tfe^-i^s {(p„ - 9sA«t a ) 2 - e f A»(p„ - g s Ap a ) 

(11) 

Eq.(ll) shows that the chromoelectric EDM term disappears, and instead 
the coefficient of the electric EDM term now changes into a new form. Fur- 
ther, the last term which is similar to the QED case of eq.(6) is generated. 
This term gives rise to the neutron EDM in the second order perturbation 
theory in the same way as eq.(7). 

Here, it should be noted that the two representations of the EDM operators 
are indeed equivalent when we estimate their expectation values with exact 
wave functions. However, when we consider the EDM of any composite 
systems like neutrons, the difference between the two operators becomes 
important since we cannot normally obtain the exact wave functions of 
neutrons. Employing the new EDM operators, we can calculate the EDM 
of neutrons in a transparent and economical fashion since it can be evaluated 
in the first order perturbation theory. 
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4. Neutron EDM 

Now, we want to estimate the neutron EDM based on the lagrangian density 
obtained in the previous section. The EDM hamiltonian which corresponds 
to the new EDM lagrangian can be written as 

H edm = - E { (d f ~ g^V) 70S • E - % 0l5 (pi - 2g s A a • pt a + g 2 s (Ap a ) 2 ) 

(12) 

where we dropped the term which is proportional to the electromagnetic 
vector potential since the neutron state does not contain any photon 
states in the ground state. Also, the magnetic term is not included in 
eq.(12). 

Now, we estimate the first term of eq.(12) in the first order perturbation. 
Using the SU(6) quark model, we obtain for the neutron EDM, 



\{d d - ^dd) - \{d u - ^-dv 
3 2g s 3 2g s 



N R (13) 



where Nr denotes the overlapping integral, and Nr = 1 for nonrelativistic 
quark models and Nr ~ 0.5 for the MIT bag model [24]. 

Next, we evaluate the contribution to the neutron EDM from eq.(12) in 
the second order perturbation as given in eq.(7). Since there is no \ state 
nearby in the neutron excited states, there is no enhancement, contrary to 
the atomic cases [21]. 

Here, we carry out a naive estimation of eq.(7) with the MIT bag model 
[24], and therefore, we employ the closure approximation. In this case, 
eq.(7) becomes 

d^ = ^<m e E {U-^d,)^-^ 

i=d,d,u j=d,d,u \ \ " s / 

- %iH 3) {pf - 2g s A a ■ pVq + gliA^f) } \N > (14) 

where AE is an average excitation energy of neutron | states which should 
be of the order of nucleon mass M^. denotes the electric field which 
the j—th quark feels inside the bag. After some numerical calculations, we 
obtain 



3R n AE 



2e 2 d id d - 7^-d d j + el ld u - ^~d v 



4:Rq / 2 7 

— a~F [Zedm d dd + e u rr 



9g s AE 



(15) 
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where m d (m u ) denotes the mass of the d (u) quarks. Now, we take the 
AE ~ Mn and Ro ~ — with m n the pion mass, and therefore, the total 
neutron EDM can be written as 



dn 



N, 



R 



2m n 
3M, 



N 



2 e 2 d (d d - 7^d d j + e 2 u ld u - ^-d u 



9g s m n M_ 



N 



(2e d m 2 d d d + e u m 2 u d u ^ . 
(16) 



Since the second line of eq.(16) is much smaller than the first line contribu- 
tion, we can practically neglect the second order perturbation contribution 
to the neutron EDM as long as the first two terms (the first line) survive. 



5. Possible suppression of neutron EDM 

As can be seen from eq.(16), the neutron EDM must be significantly sup- 
pressed if the following relation holds for / flavor quark, 



df = —d f . 



(17) 



Recently, Chang, Keung and Pilaftsis [22] presented the two-loop calcula- 
tion to the EDM in supersymmetric theories. They show that the electric 
EDM df and the chromoelectric EDM df of a light fermion / can be given 

as 

(18a) 



3a em Rfmf sr n 2 u 

a q=t,b 
rl - a ^ R f m f V H 

f ~ 9s l2%^ Ml q 

a q=t,b 



(186) 



where a em and a s are the electromagnetic and the chromomagnetic fine 
structure constants, respectively. Also, Rf is given as Rf = cot/5, and M a 
denotes the tree level mass of the CP-odd Higgs. H q is defined as 



\ Ml 



f(^ 

\ M„ 2 



(19) 



where F(z) denotes a two loop function as given in ref.[22], 



F(z) = f 
Jo 



dx 



x(l — x) 
z — x{l — x] 



In 



x(l — x) 
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Also, £ ? 's are the CP violating couplings which are given in the MSSM 
scheme [7]. 

Now, we want to see the relation between df and df. From eqs.(18), we 
find 

Since the value of ( — \ s c l ose to unity at the electroweak scale, the 
V a s J 

parameters in H q will control the magnitude of the curly bracket of eq.(20). 
If the value of the curly bracket is unity, then the relation of eq.(17) holds. 

At the present stage, we do not know whether the parameters that appear 
in eqs.(18) can be estimated reliably or not in the supersymmetric theories. 
If these parameters can be determined reliably, it would be very interesting 
to examine the validity of the relation. 
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